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Abstract 



We have advanced our previous static theory of polymer entanglement in- 
volving an extended Cahn-Hilliard functional, to include time-dependent dy- 
namics. We go beyond the Gaussian approximation, to the one-loop level, 
to compute the frequency dependent storage and loss moduli of the system. 
The three parameters in our theory are obtained by fitting to available ex- 
perimental data on polystyrene melts of various chain lengths. This provides 
a physical representation of the parameters in terms of the chain length of 
the system. We discuss the importance of the various terms in our energy 
functional with respect to their contribution to the viscoelastic response of 
the polymeric system. 
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In a previous paper, we developed a static field theory of polymer entanglementB, in 
which we introduced a non-local attractive term, in addition to the usual excluded volume 
term, that models resistance to the motion of polymers due to entanglement. Starting with 
this energy functional, we were able to use Renormalization Group techniques to describe 
the onset of entanglement as the average molecular weight is increased to a critical value. 
The onset of entanglement may be described as a cross-over phenomenon, characterized by 
the effective diffusion constant going to zero as the transition point is approached, thereby 
implying critical slowing down. We pointed out the existence of experimental evidence to 
support the theory. 

There have been several numerical approaches developed to understand the viscoelastic 
response of polymersil^i. It is of interest to see whether an alternative theory of viscoelas- 
ticity of polymers using continuum concepts can be developed. Our previous theory, being 
static in nature, clearly needs to be extended if one is to study the time-dependent response 
of polymeric systems. The chief purpose of this paper is to lay down the foundations of 
a time-dependent field theory of entangled polymers, primarily through a comparison with 
experimental results on the linear viscoelastic response of polymer melts. In future work, 
we shall probe the time- dependent approach to entanglement of polymeric systems, as the 
molecular weight is increased to a critical value. In this, we have in mind an analogy with 
dynamic critical phenomena, where one of the quantities of interest is the frequency de- 
pendent diffusion constant, and the manner in which it scales to zero as the transition is 
approached. We will also compare the results of our theory briefly to the results of the 
standard reptation theory approach!. 

The time-dependent internal energy functional U which extends our previous static 
theoryS can be written down in a straightforward manner in terms of an energy density 
u: 

U = J d 3 s J dt u(s,t) 
dc(s, t) 



(3u = c(s, t) 
+ 



dt 

a \ <9c(s, t) <9c(s, t) 
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where c is the number concentration of the polymer strands, and 5 2 = V2a, k is Boltz- 
mann's constant and T is the temperature. This model is an extension of the standard 
Cahn-Hilliard approachi. The first two terms in Eqn. ([]]), taken by themselves, are easily 
seen to lead to the time-dependent linear diffusion equation, with a plays a role analogous 
to that of a diffusion constant. The third term in the equation represents the standard 
excluded volume interaction. We will discuss a in greater detail shortly. The final non-local 
attractive term represents the fact that when polymers become entangled, there will be in 
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general a resistance to their movement. Similar approaches can be found in the literature, 
where authors have taken activation energies to represent entanglement!. The form of our 
interaction term is novel, however, and we have discussed in the previous paper! the ra- 
tionale for choosing the parameter a in the model is as given in that paper, since it leads 
to a number preserving Euler-Lagrange equation. It is worth noting that even standard 
treatments involving the excluded volume term are generally restricted to the static case, 
whereas we have considered here an extension to the frequency-dependent case as well. 

The time and space co-ordinates in the equation are dimensionless, the scales we have 
chosen being lo~ x and A respectively. The dimensionless parameter a can be written as 
D/[\ 2 oj c ], where D is the diffusion constant which scales as the inverse square of the molec- 
ular weight in the entangled state. Since both A and u c could possess a molecular weight 
dependence of their own, it follows that a can display a dependence on molecular weight 
different than the scaling behavior of D. 

Since we wish to study polymer melts undergoing shear experiments, where they are in 
contact with an energy reservoir (at constant temperature), the energy which is conserved 
is the Helmholtz free energy S being the entropy and T the temperature. 

The entropy is given approximately as!: 

S= J d 3 s J dts(s,t) 
s(s, t)m c(s, t) ln[c(s, t)] 

/2 /3 /4 

(1 + c') ln(l + c')« c' + — - — + — (2) 
v ; v ; 2 6 12 v ' 

The last of these equations indicates an expansion around A -3 , which is identical to one 
in the units we have chosen. The linear terms will be ignored following convention, as they 
can be absorbed into the chemical potential /i, required for number conservation. In the 
mean field approximation, /i = 0. In what follows, we shall drop the primes on the number 
concentration variable. 

The goal of this communication is to compute the linear viscoelastic response of a poly- 
meric system. This can be done following closely the analysis in our previous papern, to 
obtain an expression for the frequency dependent stress a{uj): 

a(uf)= —iuj CS{k = 0,u)e(uj) 

where S(k,u) is the two-point correlation function for the system. Note that since we 
chose to take temporal Fourier transforms with respect to exp(— iut), our sign convention in 
the first of Eqns.(|3|) is opposite that in standard literature!. In general S(k,u) is given by: 

-l 



S(k,u)= (S \k,u) - E(fc, 

S (k,uj)= (—iuj + SQ 1 (k) S j 

s (k)= (l + V2ak 2 + 2a 2 k 2 /(l + v 7 ^ 2 /"))^ ~ (l + ak 2 

a= 2V2a (4) 
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where, as usual £ denotes the self-energy. From Eqn.(f|), we see that a plays the role 
of a diffusion constant. But bearing in mind the discussion below Eqn.Q, we expect its 
dependence on molecular weight to be different than the conventional diffusion constant D, 
due to the manner in which we have scaled our variables. 

From Eqns.(|^) and (f|), we see that we need to evaluate the correlation function in the 
long wavelength limit. We can perform this calculation using a perturbation expansion 
with respect to the nonlinear terms, using standard methods from field theory. These are 
elementary extensions of the methods detailed in Ref. 1. The vertices we obtain from 
Eqn.(Q) are depicted in Fig.llR The diagrams which we consider are depicted in Figures 2 
and 3. It is easy to show that the tadpole and bubble diagrams vanish identically and the 
only surviving lowest order diagram is the setting-sun diagram, whose contribution can be 
shown analytically to be: 



S 2fo (u;, k = Q,w m )= - J— y (2^)3 S o(k 7 uj') S (k,u - uj') 



\Jiuj — I 
327T 2 a 3 / 2 



In 



{u/2 + u m ) {3u/2-uo m ) 



(ui/2 — u m ) (3u/2 + u n 



(5) 



Here, u is the frequency at which the system is being sheared. In order to obtain the final 
expression for the self-energy in the long wavelength limit, we first performed the frequency 
integral. This integral is logarithmically divergent, rendered finite through the use of a 
high-frequency cut-off u m , using the conventional procedure from field theory. u m may be 
interpreted in the usual field-theoretic sense as an upper frequency limit below which the 
continuum theory is not valid. The dependence of the final answer on the cutoff u m is 
logarithmic. Next, to perform the k-integration, we used the method of contour integration, 
taking care to distort the contour to avoid the branch cut implied by the logarithmic behavior 
of the integrand in k-space. Once this is done, it is straightforward to use the method of 
residues to get the expression given above. It is easy to verify that the integrand in k-space 
is indeed suitably convergent, and no additional cutoffs are necessary. 

With Eqns.(3)-(5) in hand, we fitted experimental data on polystyrene melts character- 
ized by a very low polydispersity, and a fairly wide range of molecular weights, ranging from 
just under 10 4 to about 580, 0000. The results are encapsulated in Fig. 4. Three sets of 
data are shown in the figure, chosen to represent low, medium and high molecular weight 
samples. We have performed the fitting procedure for all the sets of data provided by Onogi 
et altil. We have chosen to display three of the fits as being representative of the procedure. 
We see that the fits, while good, show signs of deteriorating slightly as the molecular weight 
decreases. The plateaus indicate the rubbery phase of the system. Our previous paper ap- 
plies in this regionS. For higher frequencies (short times), the stress is much higher, which 
could be interpreted in terms of inertia as the system starts to be strained. For frequencies 
below the plateaus (long-time behavior), one might say that the polymers eventually begin 
to disentangle, causing the stress to start decreasing precipitously. 

We began with four parameters in the theory, viz., A, a, u m and u c . In performing the 
fitting, due to the manner in which the parameters appear in the expression for the complex 
modulus, it was more convenient to use the reduced set of three parameters C, u m , a, where 
C and a have both been defined earlier. While it is true that the parameters in our theory 
had to be chosen to get the best fit with data, we find it impressive that the function given 
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in Eqn.fl^) is such that it provides the correct form for the storage modulus. In this manner, 
our theory has captured the essential aspects of the linear viscoelasticity of polymer melts. 
From the values obtained for the three parameters, we were able to perform a least-squares 
fit, yielding the following representations as functions of the molecular weight M n . 

3.1 x KT 13 Ml 1 (dyne -s- cm- 2 ) 
b M% (s- 1 ); b = 0.11; 61 = -1.53 

d M^ 1 ; d = 0.21 x 10" 13 ; d x = 1.35 (6) 

Note the scaling forms for these representations. First of all, as will be shown shortly, 
C can be identified with the static viscosity, and we found the scaling exponent 3.7 to be 
reasonably close to the value of 3.4 given in the literature. To obtain this representation, it 
was useful to plot the fitted values of C against the molecular weight M N on a log-log plot. 
The cut-off frequency u m scales inversely as approximately the 3/2 power of the molecular 
weight. Finally, as advertised earlier, the dependence of a on M n is different than that of 
the true diffusion constant D. The reason for this is that a is proportional to DX~ 2 u~ , so 
that the product \ 2 u c , which has the units of a diffusion constant, goes approximately as 
~ M" 10 / 3 . It is further easy to show that A ~ M~ 2//5 , and tu c ~ M^ 2 . 

As one might perhaps expect, the entropy terms in our energy functional have a dominant 
effect on determining the linear viscoelastic behavior of polystyrene melts. The nonlocal 
attractive term, which models the resistance to the motion of entangled polymers has a less 
pronounced effect on linear viscoelasticity. This is consistent with our earlier calculations in 
the static regimeB, where we found that the nonlocal attractive term has a more profound 
effect on determining the renormalized diffusion constant than the elastic moduli. We will 
tackle the frequecny dependence of the renormalized diffusion constant in future work. 

From Fig. 4 we see that the fits to data using our theory are quite good, as are those 
using the standard reptation theory!. The advantange of our theory is that while the 
reptation theory is restricted to the regime of highly entangled systems, and uses mean- 
field concepts, we can compute the effect of fluctuations using Feynman diagrams. With 
these parameteric representations, we were also able to compute the loss moduli for the 
samples. The curves gave a reasonable but only an average fit for the various samples, 
In other words, the loss moduli obtained through our procedure was not very sensitive to 
the parameters we obtained. Nevertheless, we note that this is an improvement over the 
standard reptation model approaches, which give a null loss modulus!. We also went a bit 
further, and attempted to use the above representations to calculate the storage moduli of 
another polymer melt, and found reasonable agreement, suggesting that our approach has a 
universal flavor to it. On the other hand, the results for solid polymers were quite abysmal, 
indicating that our fitting procedure must be redone for non-melts. We have been unable 
to locate data for solid polymers or other polymer melts having the same wide scope as the 
results of Onogi et al0 for polystyrene melts. 

It is easy to show that in the zero frequency limit, the viscosity is given by: 

r)o = Lim^o uj~ Im \—iujCS(k = 0, to)] = C (7) 

One similarly obtains an expression for the shear modulus in the zero frequency limit by 
considering the storage modulus. For sufficiently large values of the frequency, the experi- 
mental data show an approximately linear behavior, as does the present theory. 
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As mentioned at the beginning, the purpose of this note is to lay down the foundations of 
a time- dependent field-theoretic approach to the viscoelastic response of polymeric systems. 
Our eventual goal is to understand the dynamic approach to entanglement, analogous to 
treatments of dynamic critical phenomenatj. 

I would like to acknowledge a useful comment by Sanat Kumar concerning the calcula- 
tions. This research is supported by the Department of Energy contract W-7405-ENG-36, 
under the aegis of the Los Alamos National Laboratory LDRD polymer aging CD program. 
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FIGURES 



FIG. 1. (a) is a pictorial representation of the cubic term in A. Each leg corresponds to a 
factor of c, the field. The intersection of the three legs symbolizes a factor of 7 = 1/6, the coupling 
constant, (b) is a pictorial representation of the quartic term in A. A factor of —1/12 is to be 
inserted at the intersection. 

FIG. 2. (a) represents the tadpole diagram which is crucial in our calculations, (b) represents 
the setting sun diagram. Both (a) and (b) are second order contributions to the correlation function 
coming from the cubic interaction term, the first order corrections being null. 

FIG. 3. This figure represents 1-loop (bubble) contribution from the quartic interaction term 
in A. 

FIG. 4. This plot shows a comparison of our theory (solid line) with the experimental data 
by Onogi et al on polystyrene melts. The molecular weights are placed alongside the different 
sets. The comparison is fairly good, and appears to deteriorate slightly as the molecular weight 
decreases, indicating that the theory works better for high molecular weight melts. The plateaux 
indicate the rubbery regime for each sample 
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Chitanvis, Fig. [I] 
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Chitanvis, Fig. |3] 
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Chitanvis, Fig. f| 
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